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Abstract 

Galileon theories are of considerable interest since they allow for stable violations 
of the null energy condition. Since such violations could have occurred during a high- 
energy regime in the history of our universe, we are motivated to study supersymmet- 
ric extensions of these theories. This is carried out in this paper, where we construct 



> 

generic classes of M = 1 supersymmetric Galileon Lagrangians. They are shown to 
admit non-equivalent stress-energy tensors and, hence, vacua manifesting differing con- 

O ■ 

ditions for violating the null energy condition. The temporal and spatial fluctuations 
of all component fields of the supermultiplet are analyzed and shown to be stable on 
a large number of such backgrounds. In the process, we uncover a surprising connec- 
tion between conformal Galileon and ghost condensate theories, allowing for a deeper 
^ \ understanding of both types of theories. 

1 Introduction and Overview 

Matter in the universe is typically assumed to satisfy the Null Energy Condition 
(NEC) [1]. This is because standard two-derivative theories generically lead to the 
appearance of ghosts or gradient instabilities on NEC- violating backgrounds [2]. More- 
over, higher-derivative theories — generically associated with equations of motion that 
are of third- and higher-order in derivatives — also lead to the appearance of ghosts 
and are, therefore, catastrophically unstable [3] — [5]. 

However, in recent years it has become clear that these theoretical limitations are 
neither necessary nor, perhaps, desirable. Indeed, there are a number of cosmologi- 
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cal situations in which violations of the NEC become inevitable, including inflation- 
ary models with extra dimensions [6, 7], string gas cosmological scenarios [8] — [fl] 
(see [12] for a review), and pre-big bang [13] and ekpyrotic theories with cosmic 
bounces [14] — [50] in which the universe reverts from contraction to expansion (see [51, 
52] for reviews of ekpyrotic/cyclic theories). Furthermore, there is an important caveat 
to the theorem of [2] — namely, the existence of stable, NEC-violating higher-derivative 
theories that nevertheless lead to equations of motion with at most two derivatives act- 
ing on any field. Two classes of such theories have been studied in the literature so 
far. The first is ghost condensate theories [53], in which the Lagrangian is taken to 
be an analytic function of a scalar field <p and X = —^(d(f)) 2 . The second example 
is provided by the Galileons [54] — [64], in which higher-derivative terms are combined 
precisely so that the equations of motion have at most two derivatives acting on each 
field. Both types of theories have the remarkable property that they allow for stable vi- 
olations of the NEC and, hence, both types of theories can be used as effective theories 
to model novel cosmological scenarios. Indeed, ghost condensates have been used to 
violate the NEC [65] and, hence, enable cosmic bounce from contraction to expansion 
in the New Ekpyrotic Scenarios in [34, 35, 36]. Similarly, galileons have been used to 
devise cosmological scenarios in which the universe expands from asymptotically flat 
initial conditions [66, 67]. 

The relationship of these theories to string theory is not yet entirely clear, although 
it is interesting that (in a certain small field limit) the Galileon theories describe the 
fluctuations of a brane embedded in a higher-dimensional space-time [69] (see also [70] 
for a derivation of the Galileons via compactification of Lovelock gravity). Such sce- 
narios arise naturally in heterotic M-theory [71] — [74], for example, where five-branes 
wrapped on holomorphic two-cycles [75] — [79] can exist in the five-dimensional bulk 
space. The visible sector of such theories can contain exactly the supersymmetric 
standard model [80] — [86] and, hence, present realistic vacua to explore Galileon cos- 
mology. If ghost condensate or Galileon theories turn out to be relevant in modeling 
the dynamics of the universe in the high-energy regime, then it would seem necessary 
to consider these theories in a supersymmetric context. Quite independently, it is of 
theoretical interest to have a model that allows one to study the interplay between 
supersymmetry and NEC violation. In previous work [87], we super symmetrized the 
ghost condensate models. In this paper, we further this viewpoint by constructing 
generic N = 1 supersymmetric extensions of Galileon theories. 

In [87], we studied the M = 1 supersymmetric extension of ghost condensate theo- 
ries using chiral multiplets. As reviewed in Sec. 2, the extra scalar and auxiliary fields 
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required by supersymmetry are well-behaved for such models, while the fermionic mem- 
ber of the supermultiplet is not. Specifically, the fermion kinetic term on the ghost 
condensate background violates Lorentz-covariance — the spatial-derivative part has 
the wrong sign while the time-derivative part has the correct one. In this paper, we 
show that a manifestly super symmetric interaction can be added to this theory which 
has the property of restoring the fermion kinetic term to its canonical form. This will 
be the subject of Sec. 3. When we examine the effect of adding this term on the bosonic 
part of the theory, we find a surprise: the resulting scalar field theory is precisely the 
second- and third-order conformal Galileon theory! Hence, in rendering the fluctua- 
tions around the ghost condensate background canonical, we re-discover the second- 
and third-order conformal Galileon model. Moreover, using a field redefinition, the 
ghost condensate background is easily seen to be equivalent to the "self-accelerating" 
de Sitter solution of the Galileon theory. The real difference between these two theories 
lies in the form of the spatial gradient terms. Generally speaking, the Galileon theories 
are much better behaved with regard to spatial gradients than their ghost condensate 
counterparts, as will be discussed in Sec. 4. 

When we consider the fourth-order conformal Galileon Lagrangian, we find that 
there are now many choices in how to construct a supersymmetric generalization. Since 
these choices become vastly more numerous for the fifth- (and highest) order Galileon 
theory, we only present terms up to fourth-order in this paper. There are two reasons 
for this proliferation. The first is that, using integration by parts, one may rewrite 
a given action into one that is equivalent up to total derivatives. If we now discard 
the total derivatives and supersymmetrize the new action, we generically end up with 
inequivalent results. This is an important ambiguity regarding higher-derivative the- 
ories — that is, different theories, though related by integrations by parts (on a flat 
background), lead to different stress-energy tensors and, hence, different conditions 
for NEC violation. We provide a detailed treatment of these issues in Sec. 5. The 
second reason for this proliferation of choices is intrinsically supersymmetric. As in 
our previous paper on ghost condensate theories, we construct supersymmetric exten- 
sions of many-field terms by using a number of smaller building blocks. And for some 
terms, there are several inequivalent ways of subdividing them into separate building 
blocks. It is interesting to note that the resulting supersymmetric extensions can be 
quite different, and can now contain non-canonical fluctuations of the fermion field for 
example. This will be discussed in detail in Sec. 6. 

In Sec. 7, we discuss our results and speculate on future applications. Two appen- 
dices are included. The first provides useful formulae regarding the supersymmetric 

3 



building blocks that we are using, and the second discusses in more detail those su- 
persymmetric extensions that contain non-covariant kinetic terms for the fermionic 
fields. 



2 Supersymmetric Ghost Condensate 
2.1 A Review of Ghost Condensation 

The simplest form of a "ghost condensate" [53] arises within the context of a single 
real scalar field 4> in four dimensions. Assuming space-time is flat and non-dynamical, 
the evolution of is governed by a higher-derivative Lagrangian of the form 

C = P(X) , (1) 

where P(X) is an arbitrary function that is analytic around zero in 

The mass scale m is introduced to render X dimensionless. To simplify notation, we 
set m = 1 in most of the paper. For purely time-dependent solutions, the associated 
equation of motion is given by 

d 



dt 

Clearly, <p = const, is a solution. However, (3) also allows for solutions with arbitrary 
constant X, that is, 

4> = ct, (4) 

where c is a constant. Although in this paper space-time is taken to be non-dynamical, 
we note that in a cosmological context the equation of motion on a Friedmann- Robertson- 
Walker background becomes 

^ (a 3 P x <]>) = , (5) 

where a(t) is the scale factor of the universe. For a generic choice of P(X), this implies 
that <j) must redshift as the universe expands. However, there is one key exception: if 
P(X) has an extremum at some X = c 2 /2, then (f) = ct is a solution to (5) indepen- 
dent of the behavior of a(t). Moreover, this solution is an attractor on an expanding 
background — small departures away from the extremum are driven to zero by Hubble 
friction. This solution, cf) = ct, spontaneously breaks Lorentz invariance and is called 
a ghost condensate. 
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Returning to a flat background and expanding in fluctuations 

cf) = ct + S(f)(t,x) (6) 
around a ghost condensate, to quadratic order in 5(f) the Lagrangian becomes 

£quad = XP XX • («ty) 2 - • S&'Sfo . (7) 

As a result of Lorentz-breaking, the coefficients in front of the time and spatial deriva- 
tive terms are unequal. We see that the condition for the absence of a ghost is 

XP xx >0, (8) 

which is automatically satisfied close to a local minimum of P(X). (For a general X = 
const, solution, the ghost-free condition is XP XX + P >x /2 > [88].) We, henceforth, 
assume this is the case. However, the vanishing of the second term in (7) is troubling, 
since it clearly signals that the ghost condensate is on the verge of a gradient instability. 
Can this potential instability be removed? Happily, the answer is affirmative, although 
it requires introducing higher-derivative terms, such as [53] 

(9) 

into the Lagrangian that are not of the P (X) type. Such corrections are expected from 
an effective field theory point of view. Because this term involves two derivatives per 
field, the background <fi = ct clearly remains a solution. However, (9) does affect the 
gradient term of the fluctuations, giving rise to the dispersion relation uj 2 ~ k 4 /M 2 . 
For large enough mass M, this higher-derivative term can be consistently treated as a 
small correction. Be this as it may, the question of temporal ghosts and/or gradient 
instabilities in ghost condensate theories is an important one, and will become even 
more important in the supersymmetric context. 

Before proceeding, we introduce the following simplification. Sufficiently close to the 
ghost condensate point, P(X) is approximately quadratic. Without loss of generality, 
one can rescale the field (f> so that the minimum lies at X = 1/2 (corresponding to 
c = 1) and write the prototypical ghost condensate action as 

£ = -X+X 2 = +i(d0) 2 + i(d0) 4 . (10) 
The quadratic Lagrangian (7) now becomes 

£quad = m 2 ~ • SPB^i . (11) 

We will use Lagrangian (10), which contains all of the essential physics, to supersym- 
metrize ghost condensate theories. 
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2.2 Supersymmetric Ghost Condensate 

In [87], we presented an N = 1 supersymmetric extension of the bosonic ghost con- 
densate theory in (10). To do this, consider a chiral superfield 

$ = a + iOa^oA^ + -eeeeuA + oof + V29ip - -^=ee^a"e , (12) 

'4 y/2 

with the complex scalar A(x), the auxiliary field F(x) and the spinor if) a (x) being 
functions of the ordinary space-time coordinates x^. Spinor indices which we do not 
write out explicitly are understood to be summed according to the convention ip9 = 
ifj a a and ~4)9 = tpaG a - The complex scalar is chosen so that 

A = -L(cf> + i x ) , (13) 

where <fr is the real field of the bosonic condensate theory. The imaginary component 
X is a new real scalar degree of freedom, introduced into the condensate theory by 
super symmetry. That is, 4> is taken to be the lowest component of the M = 1 chiral 
supermultiplet (cf), X ,if>, F). 

It was shown in [87] that a supersymmetric extension of the prototypical ghost 
condensate Lagrangian (10) is given by 1 



(14) 



where \ ge gQ indicates taking the #6^f9-component of a superfield. (Here and throughout 
the paper, derivatives are understood as acting only on the nearest superfield, unless 
noted otherwise. For example, D<S>D<S>D&D& = (D<$>)(D<S>)(D&)(D&). Similarly 
for space-time derivatives acting on component fields.) In terms of component fields, 
(14) becomes 

- ^^ )V l -(^o^-^V) + ... (15) 

where we display terms to quadratic order only in ip and set F = 0. Note that for 
X = if) = 0, this expression exactly reduces to Lagrangian (10). It is in this sense that 
£SUSY ^ supersymmetric extension of the prototype bosonic condensate theory. 



1 As discussed in [87], there exists a second (inequivalent) supersymmetric extension of X 2 , which however 
leads to the exact same issues with the fermionic fluctuations as those discussed below. Our discussion, and 
the cure proposed below, are thus general. 
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Since x always appears at least to quadratic order, it can consistently be set to zero. 
Thus the equations of motion can be solved by the ghost condensate 



ct, x = 0. 



(16) 



The classical fermion solution is, of course, zero. That is, the Lorentz-violating ghost 
condensate continues to exist as a vacuum of the supersymmetrized theory. 
Setting c = 1 and expanding in fluctuations 



4> = t + 5<t>(t,x), x = Sx(t,x), ift = 6ift(t,x) 
around this vacuum, we find to quadratic order that 



(17) 



rSUSY 
^quad 



o- 



+ 



- (6ip fi o°6i/> - 5ipa° <fy i0 ) - - (Sip^Sip - Sift^Sift^ . (18) 



The first line reproduces the standard result (11) for the single (ft field ghost condensate, 
as it must. That is, the time derivative term is ghost-free but, at the minimum of P(X), 
the spatial gradient term for 5 (ft vanishes. As discussed above, higher-derivative terms 
of the form (9) cure the potential gradient instability in the bosonic theory and stabilize 
the dispersion relation. Can one find a supersymmetric generalization of these terms? 
In [87] we showed that this can indeed be done. The simplest such example is 

— \rD®D<S>D&D& ({D,D}{D,D}(<S> + &)) 2 

- '• V ' aaaa , quad 



-W) 2 , (19) 



where we have evaluated this up to quadratic order in fluctuations around a ghost 
condensate background. To this order, (19) does not contain Xi ift or the auxiliary field 
F at all. 

Now consider the second line in . This is the kinetic term for the scalar 

fluctuation 8x and, hence, is new to the supersymmetric theory. Note that this suffers 
from two serious problems. The first is that the temporal derivative term vanishes 
and, hence, this field is marginally a ghost. Secondly, the spatial gradient term has 
the wrong sign. Fortunately, it was shown in [87] that supersymmetric terms can be 
added to (14) that solve both problems. These are, for example, 

^D$D$M t M t ({£>, D}(& - &){D, D}(& - $)) 

— %rD<S>D<S>D& ({D,D}($ + &){D, D}(3> - &)) 2 
16 v > J e6»ee,quad 

= -2(d(ft)\dxf- {d<ft)\d(ft-dxf. (20) 
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Adding these to Lagrangian (14), and expanding to quadratic order around the ghost 
condensate, changes both the time and spatial gradients of x m (18) to the Lorentz- 
covariant expression 



-SUSY 
-'quad 



= . . . + {5 x f - 5x l 5x,i + ■■■ 



(21) 



This renders the x fluctuations stable, without adversely affecting anything else. In 
particular, since (20) vanishes when X is set to zero, the sum of (14) and the superfield 
expression in (20) remains a super symmetric generalization of the P{X) bosonic theory. 
It will be helpful in the next section if we analyze this result in more detail. First, note 
that adding (20) to the second term in (14) gives 



-^D<S>D<S>D&D& ({D,D}($ - &){D,D}(& - $)) 

I 2 1 

j^D$D$D&D& ({D,D}(3> + &){D,D}($ - + —D$D$D&D$ 



-2(d(f>)\dx) 2 -(d(f>)Hd<l>-dxf 



+ 



--( d( f>)\dx) 2 + (dcf>.dx) 2 



(22) 



where we have not shown irrelevant pure <f> terms or terms involving fermions. When 
evaluated around the ghost condensate vacuum (16) with c = 1, (22) reduces to 



-2(dxf - {xf 



+ 



(d X ? + (X? 



-\{dx? 



(23) 



That is, adding (20) to the second term of (14) exactly cancels the Lorentz-violating 
term. In addition, the signs are such that the resulting Lorentz-covariant kinetic term 
for x is ghost free with correct sign spatial gradient. Second, adding this to the first 
term in (14) produces the canonical normalization while leaving the correct sign un- 
changed. That is, 



i(9x)' 



+ 



-{d x f 



(24) 



which gives (21) precisely. 

Finally, consider the kinetic term for the fermion fluctuation Sip. This is given in 
the third line of (18) and, as with <5%, is new to the supersymmetric theory. We see 
from (18) that, although the magnitudes of the coefficients of the two dip terms are 
equal, the time-derivative term is ghost-free while the spatial gradient term has the 
wrong sign. Note that this is not the same kind of gradient instability as occurs for <j>. 
There, the coefficient of the spatial derivative term is zero or small and, hence, higher- 
derivative terms can play a role in guaranteeing stability over an extended time period. 
For ip, on the other hand, the coefficient of the wrong-sign spatial gradient term is not 
small. It follows that the inclusion of higher-derivative terms, such as those in (19), is 
necessarily irrelevant. The situation for the fermion, therefore, is more akin to that of 



the second scalar whose deep wrong-sign spatial gradient had to be corrected by the 
addition of a new second order term — the sum of the two kinetic spatial gradients 
having the correct sign. However, within the context of the super symmetric extension 
of the pure P(X) theory, we are unable to find a fermionic analog of this mechanism. 
That is, the fermion kinetic spatial gradient term has the wrong sign! 

As discussed in [87], it is unclear whether or not this is physically unacceptable. 
This will be explored elsewhere [89]. In this paper, we ask a different question: by 
modifying the bosonic theory so that it is no longer purely a P{X) theory, can one find 
a supersymmetric extension that is free of both ghost-like and gradient-like instabilities 
in all of its component fields? The answer, as we will see, is yes, and leads to another 
interesting class of higher-derivative Lagrangians — the conformal Galileon theories. 



3 Curing the Fermion Gradient Instability 

To solve the gradient instability problem for the fermion, we proceed by analogy with 
the x scalar. That is, 1) we find a supersymmetric interaction which, when added to 
the second term in (14), cancels the Lorentz-violating part of its fermion quadratic 
terms — rendering the fermion kinetic term Lorentz-covariant with the correct sign — 
and 2) we add this to the first term in (14) to canonically normalize the coefficient. 
However, there is one important caveat. As stated above, our attempts to do this with 
precisely the two terms in (14) failed. To solve this problem, it turns out that one must 
make a mild modification of each of these terms — a modification that, however, does 
not reduce to the pure P(X) theory, or even the generalized P(X, 4>) theory discussed 
in [87]. 

With this in mind, recall from (14) and (15) that 



16 



= I(^) 4 -l(^) 2 (t^-V^) 

- + - (25) 



Here and henceforth in this section, we drop irrelevant terms containing x an d set 
F = 0. Let us now modify this term to 



4($ + $t)4 



■^MAV^ - ^V) + ■■■ (26) 



To the order we are working, the only effect of (<& + 4 is to multiply expression 
(25) by an overall factor of 0~ 4 . Furthermore, setting tp = reduces (26) to X 2 /4> 4 . In 
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other words, this modified term is a super symmetric extension of the P(X, 4>) theories 
discussed in [87]. When evaluated on a ghost condensate background, the first fermionic 
term in (26) remains Lorentz-covariant, while the last term explicitly breaks Lorentz 
invar iance. 

Can one find a supersymetric interaction that will exactly cancel this Lorentz- 
violating fermion kinetic term? Consider 
-1 



24($ + $t); 



[D§D>$>D 



! $t + h.c.) 



60 3 



1 

i 



□ 



120-' 



where we work to quadratic order in the tjj fluctuations. (Useful intermediate steps in 
evaluating the above expression can be found in Appendix A). An important technical 
fact is that, while the first four terms are contained in the component expansion of 
(D<f>D<f>D 2 & + h.c .) \qqqq, the last term arises due to a contribution from the prefactor. 
This did not occur in (26) which, to the order that concerns us, was simply multiplied 
by a factor of 4>~ A . Here, however, the prefactor is significant and must be included to 
solve the fermion gradient instability problem. Integrating the second and fourth terms 
by parts, and dropping all interactions that vanish on a ghost condensate background, 
we find that (27) dramatically simplifies to 
-1 



24($ + $t)3 



(p$D$D 2 & + h.c.) 



1 



+ 



- */'»"*'•") ■ (28) 

Note that the fermion term is simply —1 times the Lorentz- violating last term of (26) - 
a fact requiring, amongst other things, the <&+<!> t prefactors in both (26) and (28). Also, 
when setting the fermion to zero (28) reduces to — □ 0(<90) 2 /60 3 , which is manifestly 
not of the P(X, <p) form. Instead, we recognize this as the cubic term of Galileon 
theories (more precisely, conformal Galileon theories, as we will see shortly). The fact 
that (28) goes beyond the P(X, (f) form is consistent with our earlier conclusion that 
the fermionic instability could not be removed within the context of supersymmetric 
ghost condensates! 

Adding (26) and (28) together, the Lorentz-violating fermion term exactly cancels 
and one obtains the Lorentz-covariant fermionic Lagrangian 



1 



24($ + $t)3 
1 



D$D§D 2 & + h.c.) + 



4($ + $t)4 



6^ 3 



-□ 



+ 



+ 



(29) 
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Integrating twice by parts, the first term can be expressed as 
It follows that 

, 24(* + *t)3 {p*D*&& + h.c.) + -^l^DtDtDtWtf 

= w) 4 + 1^2 ( w) 2 - 1^2 ( Q ^) 2 - i^i m\^i> - V^VV) +(31) 

Three fundamental conclusions can be drawn from (31): 1) the fermion kinetic term 
is Lorentz-covariant and, for any purely time-dependent background, of the correct 
sign — that is, ghost-free with correct-sign spatial gradient; 2) the first term is simply 
X 2 /30 4 and is manifestly of the P(X, (ft) type; 3) the remaining (ft terms are of a 
different differential form and not of the P(X, (ft) type. Thus, by moving away from 
purely P(X, (ft) theory we have solved the problem of the fermion gradient instability. 

As with x, one must now add this equation to the first term in Lagrangian (14). 
Since canceling the Lorentz-violating fermion kinetic term required a modification of 
the higher-derivative operators, we must also appropriately modify the first term in 
(14). Clearly, this requires multiplying the ####-component of — by 1/0 4 . Al- 
though naively one might think this would be accomplished by the expression 

_4$$t/($ + 

<3?^) 4 , the correct result is more subtle, as discussed in [90]. Defining 

*<•••*> -stiTW (32) 



-*(*. *') «*, = sir W) 2 - ^rM^"* - ■ ( 33 > 



the appropriate modification is given by 

where we suppress irrelevant % and F contributions. The first term is just —X/(ft A and 
hence of the P{X(ft) form. Although not strictly necessary, we choose to add (31) to 
(33) in such a way that the X dependent contribution for (ft takes the canonical ghost 
condensate form (10). This will be the case if one takes the complete Lagrangian to 
be (33)+3x(31): 

1 3 

^ ( ^ + ^(^ + ^( W)2 -^(D^ 



"204 ( 1 + o W) ) W.XV' - V^Vv) + . . . 



_L(_ X + x 2 ) + {(d,d u (ftf - (Ucftf) + J- (y,^ - + ... (34) 
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The first bracketed term is of the P(X, 4>) type, whereas the second group of scalar 
terms is not. Be that as it may, the ghost condensate <p = ct with c = 1 is a vacuum 
solution of the equations of motion. The coefficient of the fermion kinetic term in the 
final line of (34) has been evaluated in this vacuum. As promised, the fermion kinetic 
term is ghost free with correct-sign spatial gradients and, rescaling ip — > V2ip, has 
canonical normalization. By canonical we mean that the ratio of the fermion kinetic 
coefficient to the coefficients of the 4> and x kinetic terms is the same as in a standard 
supersymmetric Lagrangian. Note, however, that in the present case all three kinetic 
terms are multiplied by the common prefactor l/4> 4 = 1/i 4 . For completeness, we point 
out that similar vacua can be achieved by choosing the Lagrangian to be — C2 x (33)— C3 x 
3x(31) for any negative coefficients C2,C3 2 . In this ghost condensate solution 

4> = ct still exists, but with c = \Jc2 / C3 . 

Although the equations of motion continue to have a ghost condensate solution 
of the form (16), the Lagrangian (34) is not a supersymmetrized ghost condensate 
theory! Therefore, the price one pays to solve the gradient instability problem for 
the fermion is a modification of the bosonic part of the theory. Remarkably, (34), 
and its generalizations to arbitrary positive coefficients — C2 and — C3, is precisely the 
Lagrangian for a well-known class of higher-derivative models — the conformal Galileon 
theories — to which we now turn. 

4 Galileons and their Relation to P(X, (j)) Theo- 
ries 

Galileon scalar field theories were first discovered in the context of the Dvali-Gabadadze- 
Porrati (DGP) brane- world model [91, 92] and arise generically in brane-induced grav- 
ity models [93] — [100]. In a certain decoupling limit [54, 55, 101], the theory becomes 
local in four dimensions and describes a real scalar field ir (the brane-bending mode) 
with 

£DGp = _l ( ^__L (a ^O, + _L_. rS , (35) 

where T^ v is an external source, and the strong coupling scale A is related to the 
four- and five-dimensional Planck scales as A = \^6M^/Mp\. (Despite the conformal 
coupling to T^, the theory is nevertheless consistent with tests of gravity because it is 
screened in the vicinity of massive sources [54, 96, 101, 102]. See [103, 104] for reviews 
of screening mechanisms.) As a vestige of five-dimensional Poincare invariance, this 

2 The coefficients Ci are defined so as to conform with the standard notation used in the following section. 
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theory has two independent internal symmetries [55, 56], 



6f.1T 



6„TT = 



(36) 



where c and are constant. The first transformation in (36) is just a standard shift 
symmetry, whereas the second is called a Galilean symmetry. The latter protects the 
cubic interaction from being renormalized [55, 56]. Remarkably, despite its higher- 
derivative form, (35) leads to an equation of motion that is second-order in derivatives. 

In [57], Lagrangian (35) was generalized to include all possible interactions that 
are invariant under the shift and Galilean symmetries, and which lead to second-order 
equations of motion. In addition to the linear, quadratic and cubic terms in tt shown 
in (35), it was found that quartic and quintic interactions are also allowed. The most 
general "Galileon" theory is found to be a linear combination of the Lagrangians [57] 



^Gal, 2 
^Gal, 3 
^Gal,4 
^Gal, 5 



= (dirf 
= (dirf 



1 



(□vr) s 



TT^TT 



up 



7T"' 



+ -UlT1T^ V 1T 



(37) 



where we have set the associated mass scales of each term to unity to simplify notation. 
As with the cubic term, these interactions are protected by non-renormalization the- 
orems. The construction stops with £oai, 5 — no higher-order interactions can satisfy 
the simultaneous requirements of shift /Galilean invariance and second-order equations. 

The symmetries in (36) can be promoted to a subgroup of the conformal group, 
with infinitesimal dilation and special conformal transformations acting respectively as 



6 c TT = C {l+X^d ll 1T) 
6 V TT = - d^TT Q 



-V^ 2 - 

2 



r r r — (v ■ .r ).!■'' 



(38) 



In the limit of small it, these reduce to (36). The unique Lagrangians invariant under 
these symmetries and leading to second-order equations of motion are [57, 69] 

£ 2 = -^(dirf 

Cs = -i(5vr) 2 Dvr-i(9vr) 4 



-2tt 



{dirf 



-(□vr) 2 + -7r^7r 
2 y ' 2 



/(// 
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+ I( 9vr )2 D7r _ Itt ^„ _ A( 9 vr) 4 



Cn = e 



-i(Dvr) 3 -7r^ W ^ 



+ 1^^^ + ^(dvr) 2 ^) 2 - ^dirfn^n^ 



(39) 



The class of theories obtained by taking general linear combinations of these terms 
are called conformal Galileon theories. To compare Galileons to P(X, (ft) theories, it is 
useful to change variables to 

(40) 



= e 



The above Lagrangians then become 



C 2 = 
£ 3 = 

U = 



6 



703—- W 



6 ' 3 (30)* 



703 (W' 



560 4 



+ 



(41) 



where the second versions of £3 and £5 follow from integration by parts. 

Note that although the bosonic Galileon Lagrangians £2 and £3 were introduced for 
entirely different reasons, they are 'precisely of the form — derived in detail in Sec. 3 — 
required by a quadratic and cubic super symmetric theory to have a ghost condensate 
vacuum with Lorentz-covariant and canonical sign fermion kinetic energy. Specifically, 
the purely 0-dependent part of (33) and 3 x (31) are 

\2 



-£2 



(42) 
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and 



-± W>? + W) 2 - i(°^) 2 = , (43) 



40' 

respectively! Thus, as claimed in Sec. 3, Galileon theories arise naturally in general- 
ized super symmetric ghost condensate theories, independently of their original origin 
in [91, 55, 56, 57]. It is of interest, although somewhat peripheral to our main dis- 
cussion, to ask what scalar sector would emerge if we allowed the quadratic and cubic 
supersymmetrized ghost condensate theory to have a Lorentz- violating fermion kinetic 
term. This possibility is explored in detail in Appendix B. 
The most general Galileon Lagrangian is given by 

C = c 2 C 2 + c 3 £ 3 + c 4 £ 4 + c 5 £ 5 , (44) 

where the Cj coefficients are constant. Restricting to time-dependent fields only, it 
follows from (41) that 

1 6 6 

C = ^P(X) ; P{X) = c 2 X - c 3 X 2 + -c 4 X 3 - ^c 5 X 4 , (45) 

where X = (p 2 /2. In other words, for purely time- dependent backgrounds, the Galileon 
theory reduces to a P(X) theory with an overall multiplicative factor of 4>~ A . However, 
the spatial gradients are radically different, and much better behaved, than in the 
P(X, 4>) case. This fact has important consequences, which we discuss below. Nev- 
ertheless, the connection with P(X, <f>) theories considerably simplifies the analysis of 
time-dependent solutions, as we now demonstrate. 

It follows from the above discussion that the "de Sitter" solution 

7r dS = -ln(H t) (46) 

of the Galileon theory is simply the ghost-condensate solution 

= H t (47) 

of the associated P(X,<j>) theory. The constant Hq, as it is usually denoted in the 
Galileon literature, thus corresponds to the coefficient c of the ghost condensate, 

c = H . (48) 

For such solutions, where X is constant, the equation of motion derived from action 
(45) is 

2XP X -P = 0. (49) 
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In terms of the q coefficients and H , this reduces to 

c 2 - ^c 3 H% + ^c 4 F 4 - ^c 5 Hl = . (50) 

Interestingly, since the energy density p is given by 

P = ^(2XP X -P), (51) 

the equation of motion (49) implies that p = 0. That is, in order to have a solution 
with constant X, the energy density must vanish. This fact, which appeared to be 
coincidental in previous treatments of the Galileon self-accelerating solution, can now 
be seen to be a general requirement. 

What are the conditions for these X = const, solutions to be stable against small 
perturbations? As demonstrated in [57], Galileon theories have the property that they 
modify the spatial gradients of the (f> field in just such a way as to render the La- 
grangian for fluctuations 8(f>(t,x) covariant, despite the Lorentz-breaking background. 
Hence, it suffices to require positivity of the temporal kinetic term of the perturbations. 
Substituting 

<j> = H t + 5<p{t, x) (52) 

into (45) yields the quadratic Lagrangian 

1 • / 1 3 9 3 \ 

^quadratic = J^^^f (^ C2 ~ 2~ C 3 H + ~ 2°^) + ' ' ' ^ 

Stability requires the expression in brackets to be positive; that is, 

c 2 - 3c 3 # 2 + \c A H* - 3c 5 H$ > . (54) 

Note that the this inequality can easily be satisfied simultaneously with constraint (50) 
derived from the equation of motion. 



5 Violating the Null Energy Condition 

Galileons are interesting theoretically because they can violate the NEC 

T^n v > , (55) 

where is an arbitrary null vector, while having stable temporal and spatial fluctu- 
ations. This is no small feat since, under very general conditions, theories with two 
derivatives are inevitably plagued with ghost or gradient instabilities on NEC- violating 
backgrounds [2]. Galileon (and ghost condensate [65]) Lagrangians circumvent this 
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problem by having more than two derivatives and, hence, can have vacua with stable 
violations of the NEC [66]. This is particularly interesting for cosmological applica- 
tions. In a cosmological context, (55) reduces to p + V > 0, where V denotes pressure. 
Since H = — (p + V)/2, violating this inequality then allows the universe to bounce 
from a contracting to an expanding phase. 

In this section, we derive the conditions under which the conformal Galileon La- 
grangian (44) violates the NEC. Although this question has been studied in earlier 
work [66], here we point out important new ambiguities in defining the stress tensor, 
above and beyond the usual field theory ambiguities in T^y. Remarkably, even on a 
flat-space background (which is what we consider in this paper), T^ v is sensitive to how 
one defines the theory, including total divergence terms. In particular, two flat-space 
Lagrangians that differ only by integration by parts can have physically different stress 
tensors. 

Before proceeding, we point out an important distinction between Galileons and 
ghost condensate theories. Recall from Section 4 that for purely time-dependent so- 
lutions, <fi = <p(t), the conformal Galileon Lagrangian (44) reduces to the particular 
P(X,(f>) theory in (45). The full Galileon Lagrangian, of course, differs from the 
corresponding P(X, (f>) theory by spatial gradient terms, but these are irrelevant in 
computing the (j) equation of motion (50), its energy density (51) and the ghost-free 
condition (54). However, the gradient terms are important for computing the pressure. 
Setting the gradient terms to zero in the action does not commute with varying the 
action to obtain the pressure, as we will see explicitly in the examples below. 

A standard way to derive the stress tensor is by covariantizing the theory and 
varying with respect to the metric. Alternatively, entirely within the context of field 
theory on flat space, the stress tensor is derived via the Noether/Belinfante method. 
For theories that include up to two derivatives per field, that is, £(</>, d^cj), dfj,d u <f>), the 
Belinfante stress tensor is given by [105] 

- - \m$r* + * (to) "* + * (») ^ ■ 

(56) 

By construction, this is both symmetric and conserved. Moreover, in all cases we have 
checked, the Belinfante definition agrees with the covariantization method and, hence, 
gives the correct stress tensor that sources gravity. 

To compute the pressure, let us set \i = v = i and assume 4> = <j){t). In this case, 
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the second line in (56) vanishes and the pressure is given by 

where no summation is assumed in the second term. This clearly elucidates the dif- 
ference between Galileons and the corresponding P(X, <f>) theories mentioned earlier. 
For pure P(X,<f>), the second term is manifestly absent, and (57) reduces to the usual 
V = C = P(X,4>). For Galileons, however, the second term will in general contribute 
to the pressure, even on purely time-dependent backgrounds. 
First consider Li = — (<9(/>) 2 /2</> 4 - In this case, (57) gives 

where we have substituted (47) in the last step. Next, let us compute the pressure 
for £3. This is the simplest example that displays the integration- by-parts ambiguities 
alluded to earlier. We begin with the definition of £3 given by the second line in (41), 

£ lst version = ^U^f - -^(8^ . (59) 

In this case, □c/>(<9</>) 2 contributes to the second term in (57). Hence, the pressure 
differs from the corresponding P(X, <f>) result and 

9 Q Q 

-plst version _ _^_jA ±_ jA _ Ji_ 

This agrees with V3 derived in [66], since their definition of £3 was identical to (59). 
Now, instead, define £3 by the third line in (41), namely 

^3 nd ™ n = W) 4 " ^ (W) 2 + 77^(^) 2 • (61) 

Although this was obtained from (59) solely by integration by parts, the associated 
pressure is different! Indeed, 

T->2nd version 1 f«o\ 

1 3 ~~4t4' ^ ' 

which disagrees with (60). The resolution of this paradox is as follows. As men- 
tioned earlier, the Belinfante stress tensor gives the same answer as the covariantization 
method evaluated on a flat background. The point is that although (59) and (61) differ 
only by a total derivative, their covariant versions do not. Indeed, in going from (59) 
to (61) we have canceled the terms 

-L(^-na^)^. (63) 
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Although fully justified in flat space, such terms do not cancel on a curved background. 
Instead, they give rise to the non-minimal coupling 

(v^Cty - nv^) = --^R^d^d^ . (64) 

Even though this non-minimal coupling vanishes on a flat background, its variation 
does not! It is the contribution of the variation of this non-minimal term to the stress 
tensor that accounts for the discrepancy between (60) and (62). The lesson is that 
the stress tensor of Galileon theories, thanks to their higher-derivative nature, depends 
on the precise form of the theory in flat space. If two Lagrangians differ by a total 
derivative, then their stress tensors will agree provided that in the process of integrating 
by parts one only cancels terms that would also cancel on a curved background. For 
example, the stress tensor for £2 is unambiguous, but that of £3 and higher-order 
Galileon terms clearly are not. 

In the next section, we will see that in order to supersymmetrize £4, it is most 
convenient to use a new version of £4, related to the version in (41) by integration by 
parts. In order to avoid any confusion, we quote the results for the pressure for the 
conformal Galileon Lagrangians, as defined in (39), and then also for the version that 
we supersymmetrize. With the Galileon Lagrangians defined as in (39), substituting 
into (56) and setting n(t) = — In(Hat) as the background solution, the pressure is 

r 3 



_3_ 

9F n 2 



^ = "l/ (65) 



V 5 



21# 4 



8t 4 ' 

In particular, since our convention for the form of £4 and £5 differs from that of [66], 
our results for the pressure do not agree. We have checked that all of these agree with 
the covariantization method, evaluated on a flat-space background. Since p = on 
this background, the condition for violating the NEC for the full Lagrangian (44) is 
therefore 

V oc c 2 + h^Hl - ^-c A H* - ^c 5 # 6 < . (66) 

When we supersymmetrize these theories, we use the same form for £2 and £3, but 
for £4 we use instead the last line in (78), and, for the reasons described before, we do 
not supersymmetrize £5 explicitly. We then obtain 

P| USY = ^ , (67) 
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and in our case the condition for violating the NEC becomes 

P SUSY occ 2 + ^ 3J ff 2 + ^ 4 F 4 <0. 



(68) 



6 Supersymmetric Galileons 

Having discussed Galileons associated with a single real scalar field (ft, we proceed 
to supersymmetrize these theories by embedding <j) in an N = 1 chiral superfield 
= (0) Xi V 7 ) F). The procedure we follow is identical to that used in supersym- 
metrizing P(X, (ft) theories in [87] and employs formulae discussed there, such as the 
supersymmetry algebra 

{D a ,D 6 } = -2i<jZ & d li (69) 

and its immediate consequence 

{D, D}^{D, D}E = Sd^dZ (70) 

for any chiral superfields * and S. In addition, in writing the supersymmetric exten- 
sions of the Galileon Lagrangians we are making use of several new building blocks, 
whose component expressions have been written out explicitly in Appendix A. 

We find that the possible supersymmetric extensions of £2 and £3 are very limited. 
However, there are several options on how to build supersymmetric extensions of the 4>- 
dependent Lagrangian £4. When we reach £5, the choices on how to supersymmetrize 
become so numerous, and the corresponding expressions so large, that it becomes 
impractical — and not very illuminating — to write them out explicitly, though there 
is no obstacle in principle. Hence, we only consider the Galileon Lagrangians up to £4 
from this point on. 

6.1 £ 2 

The supersymmetric extension of £2 in (41) has already been discussed in Sec. 3. 
Defining 

*<•••*> -s(*TW (71) 

the complete supersymmetrized £2 action is given by 
£l USY =«■(*,**) 



- ~w {9 * ? - W [Bxf V F-F + W 1 *^ ~ ■ (T2) 

Note that this matches the corresponding expression in (41) when x = F = tp = 0. 
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6.2 L, 



In Sec. 3, we "discovered" the third-order Galileon Lagrangian by looking for a cure 
for the wrong-sign fermionic spatial gradient term obtained in supersymmetrizing the 
ordinary ghost condensate theory. We now examine £3 in more detail, including all 
fields of the chiral supermultiplet (<j>, %, F, ip). Working to quadratic order in all fields 
except for <f>, we find 



($ + 



— x-^ (d&D&D 2 ^ + h.c. 

f $t)3 V 



quad 



12 4i - - 2i 

-□0F*F + -(d0) 2 F*F + -73^,^'" " V^>,,) - -3 n^ v a v i> ~ WAv) 



2i 



6i 



(73) 



and 



($ + $t)4 



, quad 



1 W) 4 " |r W) 2 (<9x) 2 + ^ W • dxf 



4 



2i 



--(d<t>YF*F - -(d^i^a^ - V^W) - j^^AVo^ - 4»T^ U ) • (74) 

These can be combined to give a super symmetric extension of the £3 conformal Galileon 
Lagrangian 



r SVSY _ \_ 

^3, quad 8($ + $t)3 



D<S>D$D 2 & + h.c. 



4($ + $t)4 



40 
1 



60 



60 2 



_ —d»xd u xd M+ - -=n 



1 



20' 



■(<90) 2 F*F 



5? - - 2? 

+_(^) 2 ^,^> - ^>,„) - -□0(V V a ! > - V^</>) 

4i - - 2i 



(75) 



To obtain the bosonic part of (75), we have used (30) and the fact that, integrating by 
parts, 

1 



1 3 

73 D xW • d X) = 773 (d(f> ■ dxY 



jxV^,^ ~ 7^ W + 7^3 □0(5x) 2 - (76) 



Note that (75) reduces to £3 in (41) when x = 4> = F = 0, as it should. 

Finally, integrating the third and fourth fermion terms by parts and dropping any 
term that vanishes on a ghost condensate background (where <j>^ v = and, hence, 
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X = const.), -Cjj^uXj reduces to 

4 U q S uId,X=co„st = "^4 + + ^WfM " • (77) 

Remarkably, at quadratic order, the second scalar % does not contribute to £| USY on a 
constant X background. Another interesting feature of this Lagrangian is that, despite 
its higher-derivative nature, no kinetic term for F is generated. Hence, the auxiliary 
field can be eliminated as described in our previous paper [87] . 



6.3 £4 

While the supersymmetric extension of £3 derived above was relatively straightforward, 
the analogous construction for £4 is more complicated. First of all, for £4 in (41) 
the building blocks necessary to constructing supersymmetric generalizations are not 
manifest. It is useful, therefore, to use integration by parts to rewrite this fourth-order 
Lagrangian as 



5 ^ 3 V Y) Y Y Y^v 2Q ^ 4 V Y) 
1 . ,,0 ., 1 1 



^{dd)) 2 6^6' v 6 uv 

5(/> 3V v > v v ^ 2O0 4 



1 „ ,„ ,.0 , 1 I 



" -7^^(^) 2 ^ W) 2 + ^^^^'"^ - ^(^) 4 D^ . (78) 

The last expression consists of only three terms and is particularly simple. We will 
focus on this version, and construct supersymmetric extensions for each of its three 
terms. For the first term, consider 

= ^ 1 + ^ )2 {D,D}(D*D*){D,D}(D*W&) od5S . (79) 



In components, up to quadratic order in fields other than (ft and using integration by 
parts, this becomes 

A SU li Y tc rm , q uad = -^2W) 2 #W) 2 

+ i_^ W )2 (ax) 2 _ J_ DW) 2 (0x) 2 _ }_ w A ^ r X 

+ -L(d<t>) 2 dF ■ dF* + _L^(^) 2 a^(F*F) + ^»^ V F*F 
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+^d^^)a v d^U^ - * d ^nW)a»d»$<t>,„) , (80) 

which reduces to the first term in (78) when \ = ip = F = 0. Moreover, on a ghost 
condensate background, and dropping higher-derivative kinetic terms for fields other 
than (p, (80) further simplifies to 

A SU lst Y ter m , q uad,X=con S t = " ^W^fff* (00) 2 + J 2 {d<j>?dF ■ OF* 

-^mfMA^wt - i>o»$*). (si) 

There are three noteworthy features here. First, we see that the scalar x does not 
contribute at quadratic order on a constant X background. Second, note that (81) con- 
tains a kinetic term for the "auxiliary" field F. And third, on a purely time-dependent 
background, the fermion kinetic term becomes non-covariant since it contains only the 
time-derivative part. These last two issues will be discussed in detail below. 
The second term in the last line in (78) can be super symmetrized as 

^Stcrm = 128( ~| g f)2 ({A D}(* + &){D, D}(D*D*)LP& + h.C.) 



(82) 

In terms of (p and an d using integration by parts, we obtain to quadratic order in 
fields other than (j), 

2 

20 2 ^ ^ % " r ' ~^ ' 2<j)~'"" ! ' ~ 

+ ^d^dcpfx^ax + ^ Wx,^°x 
1 



4 U 2 S „d term, quad = ^'^(d^D^ + ^(d X f 



+ 
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(FU(f)) jfJi F* - {F*U4>)^F + (OF ■ d<t>)^F* 



+ ijod^fd^F^F) - ^ d ,(d<P) 2 F*F 



+(dF* ■ dcP),,F - (F<f>, v ) tli F*^ - (F*^),^" 

\_ „..o 1 

4(/> : 
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+ ^3 D W (Y^W,,) - 

+ ^d^id^) 2 ^^-^^) . (83) 

When x = V 7 = = 0) this is simply the second term in the final line of (78). On 
a ghost condensate background, dropping higher-derivative kinetic terms for -ip and 
integrating by parts, this further reduces to 

4 U 2nd term, quad, X=const = + ^(d<j>) A F*F - ^ (dtj> ■ 8F)(d^ ■ OF*) 



3i 

The expression above contains a non-Lorentz-covariant kinetic term for F on a time- 
dependent background, as well as both a covariant and a non-covariant kinetic term 
for ip. 

For the third term in the last line of (78), we have a choice of how to supersym- 
metrize since one can take either (<90) 4 or (<9</>) 2 D0 as the basic building block. This 
leads to inequivalent results. 

Choice 1 — based on (<9</>) 4 : The supersymmetric extension in this case is given by 



+ (d^f^^a^ - . (84) 



£! U 3fIterm(l), q uad = 64($ + g t )3 D*D*D&D&{D, D}{D, £>}($ + 



1 



(d^f + {d x f - 2(d<f>y(d x y + Hd<p ■ dxY - {d<t>fF*F 



+^-,{d(j)fdMa il Ui, - □Vxr'V) . (85) 



For x = "ip = F = 0, this gives the third term in the final line of (78). Because of the 
prefactor, quadratic fluctuations in x an d F vanish on a ghost condensate background 
and, dropping higher-derivative terms for ip, we are left with 

•^4^3rdterm(l),quad,X=const = "4^3 W) 4[=I <A 

-^Wf^tW " • (86) 
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Choice 2 — based on (<90) 2 D(/>: Here we have additional freedom since there is some 
ambiguity about how to write the supersymmetric version of the (d(f>) 2 factor. That is, 

H 3rd term(2), quad ~ 64($ + $t)3 ^ 4 i", U^? + V ftV, V t\W + ?> ) 



6»6>6»0,quad 



+a{D,D}^{D,D}^ 

x(id4>) 2 + a(d X f 



1 — — »7 — — 

+ ' □^^(^'V'V - <K>") , (87) 

where o is an arbitrary real number. Note that for \ = V' = F = 0, this also leads to 
the third term in the last line of (78). Specializing to a ghost condensate background, 
integrating by parts and dropping higher-derivative kinetic terms for ip, (87) becomes 

A SU 3^tc rm( 2), q uad,X=const = ^W) 4 ^ " ^4 W W " ^ + ^W)W 

+0^) ^^^ " ^V") • ( 88 ) 

The terms proportional to a have disappeared and, hence, one can choose a convenient 
value for a, such as a = 1. Note that the term containing i 71 *^ is also missing. The 
fermion appears both with a canonical and a Lorentz-breaking fluctuation term on 
a constant X background. Note that, contrary to all previous terms, £f^fd term(2) 
introduces kinetic terms for the second scalar both covariant and non-covariant. 

In general, one can use a linear combination of Choices 1 and 2 above. We therefore 
define our supersymmetric generalization of £4 to be 

/-SUSY _ /-SUSY , /-SUSY , /-, n /-SUSY , . /-SUSY (QQ\ 

*-4 — *-4, 1st term *-4, 2nd term t I 1 "J^, 3rd term(l) z -4, 3rd term(2) ' l ob V 

where b is an arbitrary real number. In components, and restricted to a ghost conden- 
sate background with quadratic fluctuation terms up to two-derivatives, this becomes 

4S,X=const = U{<t>)-^{d<t>)\d<f>-d X f + ^{d<t>)\d X f 

+ ^{dct ) ) i F*F + ^{d<t>fdF ■ OF* - ^{d4> ■ dF){d4> ■ OF*) 
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V 8 J 4 

(_i s _|_ 1 7 
± J ^(9^)VA(W " • (90) 

When x = V> = ^ = this reduces to the (^-dependent Galileon term £ 4 in (78), 
as it should. However, as it stands, the above Lagrangian still contains a number of 
troubling features. We now discuss how to address these. 

First, the scalar x appears with a non-covariant kinetic term. This, however, is not 
a serious problem since it can be dealt with in exactly the same manner as discussed 
in and below (20). That is, the non-covariant term can be eliminated by adding 

36 _ + f r „ ^, .,+,,„ ^, „+,\2 

eeee, quad 

= -^(d^(dcb-dx) 2 . (91) 



■ 2 ^$+&r D ® D * mw ^ ({A5}($+$t ){A d}($ 



Moreover, the coefficient in front of the covariant kinetic term for \ can De modified 
arbitrarily by adding a term proportional to 

—L^D^D^D^ ({D, D}(<S> - &){D, D}(& - $ 

= - j A {d<t>)\d X f (92) 

without affecting anything else. Hence, the stability of x can always be assured. 
Second, there are two troubling terms for the "auxiliary" field F, namely 

^(dc^fdF ■ dF* - -^{dc/i ■ dF){d</, ■ dF*) . (93) 

These act as kinetic terms for F which, hence, is no longer an auxiliary field. We 
note that this occurs because we are interested in a time-dependent background <fi(t). 
Around the usual zero vacuum, such terms would be higher-order interactions and 
would not trouble us unduly. A propagating F-field is not necessarily a problem. It 
would imply that both complex components of the Weyl spinor ip would now propagate 
- giving a supersymmetric multiplet with four bosonic and four fermionic physical 
degrees of freedom. However, in this paper we will follow a conservative approach and 
add the appropriate terms that restore F to its auxiliary field status. Consider the 
following supersymmetric terms evaluated for constant X, 



-^-L-^D<t>D<S>D&D&{D, D}D 2 <Z>{D, D}D 2 & 



(96»06>,quad 



^—^( d Af{dA*fdF.dF* 

^(d^dF ■ dF* , (94) 
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and 



2 10 ($ + $t)2 



D$D$D&D&\{D, D}${D, D}D 2 <5>\ 
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{dAf{dA*) 2 {dA ■ dF){dA* ■ of*) 



(A + A*) 2 



= -\(d<j)) A (d<f>-dF)(d(f>-dF*). (95) 
4> 

At quadratic order, these do not involve x or ^- Therefore, they can be added with 
suitable coefficients to £4 USY to cancel the unwanted kinetic terms for F, again without 
changing anything else. Thus, one can ensure that the auxiliary field remains truly 
auxiliary. 

Finally, consider the fermionic kinetic terms in (90). The first is covariant, and 
unproblematic. The second one is Lorentz-violating and, hence, undesirable. This 
term can be eliminated by choosing b = 6/5. With this choice, and adding in the terms 
just discussed, we find that a healthy supersymmetric extension of the fourth-order 
conformal Galileon Lagrangian is given by 



^ = ( 6 4($ + $t)2 ^ D}(D*D*){D, D}(D&D&) 



' {D, £>}($ + $ f ){Z?, D}(D$D$)D 2 & + h.c. 



128($ + $t)2 

" 5 x U (l + W D*D*D*W*1{D,D}{D,DH* + *t) 
+ 5 x 64(^ + gt)3 + h-c-H A D}*{D, D}& 

+ 28 ^ ^ 2 D$D$D$W&{D, D}D 2 <5>{D, D}D 2 & 



1 



{D,D}<S>{D,D}D 2 <S> 



(96) 



2 9 ($ + $t)2- 

In components, up to quadratic order in fields other than <fi on a constant X back- 
ground, this reduces to 



(97) 



It is encouraging to see that healthy supersymmetric extensions of the Galileon La- 
grangians exist, as demonstrated above. We would like to emphasize once more that, 

27 



as should be clear already from the discussion around (78), the super symmetric ex- 
tension of £4 given above is not unique. Hence it would be interesting to see, should 
a derivation of a supersymmetric Galileon theory be found in a more fundamental 
setting, precisely which form of the Lagrangian would arise. 



7 Discussion and Outlook 

In this paper, we have shown how to construct supersymmetric extensions of the con- 
formal Galileon theories. In doing so, we have uncovered a deep connection between 
Galileon and ghost condensate theories. That is, conformal Galileons can be seen as 
equivalent to ghost condensate models — in terms of the temporal gradients alone, 
the two theories are identical up to an overall factor of </> -4 — but with improved be- 
havior of the spatial gradients. This connection clarifies the role of both theories, and 
significantly simplifies the analysis of time-dependent solutions of the Galileon theories. 

In our analysis, we have encountered two important subtleties, one related to super- 
symmetry and the other inherent already in the bosonic Galileons. For the quadratic 
and cubic Galileon Lagrangians, the supersymmetric extensions are highly constrained 
and, around a ghost-condensate/self-accelerating-de-Sitter background, lead to covari- 
ant fluctuations for all fields. In contrast, for the quartic (and quintic) conformal 
Galileon Lagrangian there are many inequivalent ways to construct supersymmetric 
extensions. For some of these options, non-covariant fluctuations in some fields can 
arise, as well as kinetic terms for the "auxiliary" field. We have discussed these possi- 
bilities, and have provided an illustrative example of a completely healthy supersym- 
metric version of the fourth-order conformal Galileons, for which all fluctuations are 
covariant, and where the auxiliary field remains truly auxiliary. 

A second subtlety we encountered, and which is inherent in higher-derivative theo- 
ries, is that Lagrangians related using integration by parts generically lead to different 
stress-energy tensors and thus, for example, different conditions for violating the NEC. 
Keeping this subtlety in mind, let us now put all our results together and see if we can 
truly have a stable, NEC-violating solution of our supersymmetric conformal Galileon 
theory. For the Lagrangian 

£ SUSY = c 2 £| u SY + c 3 £ 3 SUSY + c 4 £i USY , (98) 



with £| USY given by (72), £| USY by (75) and £^ USY by (96), the conditions for hav- 
ing a) a ghost-condensate/self-accelerating-de-Sitter solution, b) stability, c) NEC- 
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violation and d) canonical and correct-sign fermionic fluctuations are 



c 2 - -c 3 # + -c 4 H 







(99) 


c 2 - 3c 3 Hq + ^c 4 Hq 


> 





(100) 


3 3 
c 2 + 2 C ^ H o + 2 C4H o 


< 





(101) 


3 o 18 4 

c 2 - -c 3 i? + y CaH o 


> 


o, 


(102) 



respectively. These can be satisfied simultaneously provided that 

c 2 < ^c 3 H% < 0, (103) 

with the value of c^Hq determined by (99). Hence, we have an example of a supersym- 
metric conformal Galileon theory that has a background solution which is both stable 
and can violate the NEC at the same time! We note that this is merely a proof of 
principle, and that using a super symmetric £5, or other choices for supersymmetrizing 
£4, will lead to a variety of such theories. For these, the conditions for violating the 
NEC and having canonical fermionic fluctuations will have to be re-evaluated on a 
case-by-case basis, but it seems likely to us that healthy supersymmetric models might 
also exist for which all the conditions mentioned above can be satisfied with c 2 > 0. In 
that case, even the ordinary zero vacuum would be stable. 
We anticipate a number of applications for our results: 

• Since the (non-supersymmetric) conformal Galileon theories can be derived as 
the Lagrangians describing the fluctuations of a brane in a higher-dimensional 
space-time, there does not seem to exist a fundamental obstacle to deriving the 
supersymmetric Galileons in a supergravity context. It will then be interesting to 
see precisely which version of the supersymmetric Galileons comes out naturally. 
Moreover, such a treatment would require one to extend the present work to 
local supersymmetry and the coupling to gravity, which will be of importance for 
cosmological applications. This derivation will appear elsewhere [106]. 

• As we discussed in detail in the paper, two Galileon actions, related using integra- 
tion by parts and dropping the surface terms, are physically inequivalent in that 
they lead to the same time-dependent backgrounds, but to different pressures. 
Hence, for two such theories, the conditions for violating the NEC are different. 
It may be that one theory allows for stable violations of the NEC, while the other 
does not. It will be interesting to investigate this situation further in a cosmo- 
logical context. Indeed, it means that in approaching a regime where the NEC is 
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violated, such as a cosmic bounce, spatial gradients must necessarily make their 
presence felt, and either allow or disallow entering into the NEC-violating regime. 
This should be the case regardless of how small the spatial gradients are initially! 
It will be of interest to see how this works out in a concrete model. 

Supersymmetric Galileons provide a fascinating theoretical laboratory in which to 
study the connections between higher derivatives, supersymmetry and violations of the 
NEC. This is at present largely unchartered territory, but, in this paper, we hope to 
have provided the basic tools necessary for mapping it out. 
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Appendix A: Useful Superfield Expressions 

In constructing the supersymmetric extensions of the higher-derivative scalar actions 
discussed in this paper, we make frequent use of a number of building blocks. We list 
them here for reference, our notation and conventions being those of Wess and Bagger 
[90]: 

D$D$ = 24)4) + W2F6^ + 4V2iV<r"04/i + mp2 ~ 469(dA) 2 

+2V26d{04)UA + 4>^a u eA, v - 2A^9^) 
+2y/2iee{F^a^6 - SFip^S) 

+eeee{AFoA - abf ■ oa - y>ny> + <9y> • dtp - 4)^a v 4) tU ) , (104) 

D 2 <5> = _4F + AV^itp^e - A69aA 

+4i9a^6F^ - 2y/ween%j) - OOOeaF . (105) 

Multiplying the first of these building blocks with the hermitian conjugate of the second, 
we obtain, up to quadratic order in the spinor ip, 

D$D$>D 2 $ Jf = SipipF* - 16\/26>V>FF* - IQV^iipa^BA^F* 
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-m^nA* - meF 2 F* + mee^^^F + me{dAfF* 

-32i6a^6A^FF* + 32^99a' / ^ >u A^ 
-8i9a^6%bipF*^ + !Qi9^^9F* - \U9^^oHF* 
+8^299 ( - 2#<t^Am dA * - 2a^^^a v 9FA^ v 

+i^9FF^ - i^OF^F* + 3iijj^9FF*^ 
+8\/2#fl(2^V>(&4) 2 + fl^cr^A^ 

+24'^,/^* - Otf>oAF* - tp^a^a^OA^F*^ 
+9999 (l6(d A) 2 a A* - 16D AFF* + WF*dF ■ dA - !6FdF*dA 
+8ii)^o u o x i}, x A, u - m^^a v i)^A^ + 8iipa"ip :l/ nA 
+8i^a ll U^A^ + AipDipF* - Ad^j ■ d^F* + \^ ^ o v ^ ^F* 
+4^(7^ Y^F* - 4tP^a u a^F* v - 2^DF*) . (106) 

We also make frequent use of 

($ + $t)* = ( A + A*) k + kV2(A + A*) k ~ 1 (e7p + 9ij) + k68(A + A*) k - 1 F 

+k89(A + A*f- l F* + ki0a»6~{A + A*) k ~ x {A^ - A*^) , (107) 

where we have dropped the top component as well as terms quadratic and higher in 
fields other than <p. 

Appendix B: Non Lorentz-Covariant Fermion Ki- 
netic Terms 

In Sec. 3, we showed that the exact linear combination of (28) + (26) not only gives the 
conformal third-order scalar Galileon Lagrangian, but also results in a Lorentz-invariant 
fermion kinetic term. In this Appendix, we generalize this analysis by allowing for a 
more general linear combination. The fermion kinetic term now breaks Lorentz invari- 
ance, and the resulting generalized Galileon theory is only invariant under dilations but 
not special conformal transformations. 

Instead of (31), consider the more general expression 

(28) + (1 + A) x (26) = ^ i0 tf + ^(W) 2 - ^C^) 2 + A W) < 

-^i^ m 2 (t^ - - ^MAV^ - ^V) + ■■■ 

(108) 
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where A is a constant. It follows that 3 x (31), the higher-derivative term entering 
expression (34), is now replaced by 



c A - 3?( 4^4 A W ) 2 (V>^ - Y^Yv) - ^4>AAVo^ - ^V) + ■■■ (109) 



where 



3A 

Ca^-Cs + ^W) 4 , (110) 



and £3 is the third-order conformal Galileon Lagrangian in (41). 

For A = 0, the last term in (109) vanishes, and the fermion kinetic term is Lorentz- 
covariant with the correct sign on the ghost condensate background. Furthermore, 
Ca=o = — £3 is the standard conformal Galileon action discussed in Sec. 4. As men- 
tioned in Sec. 4, £3 is invariant, up to a total derivative, under the infinitesimal dilations 
and special conformal transformations given in (38). We now show that in the more 
general case, when A / 0, only the dilation symmetry survives. 

To see this, we work in terms of tt and decompose £3 as 

C 3 = C 3A + C 3B , (in) 

with 

C 3A = ~(dir) 2 n 7T ; £ 3B = -I^vr) 4 (112) 

being of order 0(7r 3 ) and 0(7r 4 ) respectively. In this notation, the generalized La- 
grangian Ca can be written as 

Ca = ~C 3A - (1 + 3A)£ 3 b • (113) 

First consider the dilation transformation 5 c ir = c(l + x M <9 M 7r). Since each tt in (112) 
is acted on by at least one derivative, the relevant variation is 

5 c d^TT = cc^(x Q cU) . (114) 

This preserves the order in tt and, hence, the variations of £3,4 and £33 cannot cancel 
against each other. Instead they must be separately invariant (up to a total derivative) 
under this transformation. It is straightforward to check that this is indeed the case. 
And since Ca is just a linear combination of these two terms, it too is dilation invariant. 

Now consider the infinitesimal special conformal transformation, given by the sec- 
ond equation in (38). The relevant variation in this case, 

Svd^TT = V, l -d ll (d a TT (Vx 2 - (V ■ X)X^ ) , (115) 

32 



has both Oth- and lst-order contributions in it. Thus, unlike the previous transforma- 
tion, the variations of the two terms in (112) can cancel against each other, so that 
neither need be a total divergence. For example, consider 

5 V C 3B = -v^d^idir) 2 + 0(tt 4 ) . (116) 

Since 6~ v Cza is at most cubic in it, the C(vr 4 ) term in (116) must be a total derivative, 
which is indeed the case. However, the first term need not a total derivative, and, in 
fact, it is not. This is most easily checked by defining an action — / d 4 x v^{d^n){dTr) 2 . 
If the integrand were a total derivative, then the variation of this action would vanish 
identically. Instead we find the non-zero result 

-— / d 4 x^(a„7r)(a7r) 2 = 2^(5„7r)nvr + 4^(a a 7r)(a„a a 7r) / . (117) 

07T J 

It follows that Lzb is not by itself invariant under (115). The invariance of £3 relies 
on a cancellation between the first term in (116) and the 0(ir 3 ) part of the variation 
of Csa (the C(7T 2 ) term in 8 v Cza is a total derivative). 

The immediate corollary is that Ca is not invariant under the transformation (115) 
for A^O. In other words, when the fermion kinetic term is not Lorentz invariant on 
the condensate background, the purely 7r-dependent part of Lagangian is still invariant 
under dilations (114) but breaks special conformation transformations (115). This 
generalized class of "detuned" Galileon theories, and their supersymmetric extension 
discussed above, admit a ghost condensate vacuum and are potentially interesting in 
their own right, such as for cosmological applications. We will explore their properties 
elsewhere [89]. 
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